Properties of quasi two-dimensional condensates in highly anisotropic traps 



O 

o 

(N 
-i— > 
O 

O 



5h 



c 

o 
o 



> 

00 

^1- 



^1- 
o 

-I— > 



I 

c 

o 
o 



X 



G. Hechenblaikner, J. M. Krueger, and C. J. Foot 
Clarendon Laboratory, Department of Physics, University of Oxford, 
Parks Road, Oxford, 0X1 3PU, 
United Kingdom. 
(Dated: February 2, 2008) 

We theoretically investigate some of the observable properties of quasi two-dimensional conden- 
sates. Using a variational model based on a Gaussian-parabolic trial wavefunction we calculate 
chemical potential, condensate size in time-of-flight, release energy and collective excitation spec- 
trum for varying trap geometries and atom numbers and find good agreement with recent published 
experimental results. 



I. INTRODUCTION 

Bosc-Einstein condensation of dilute atomic gases has 
been achieved in a variety of magnetic and optical dipole 
force traps with different geometries. There is consider- 
able interest in studying the properties of these ultra-cold 
gases under conditions where the confinement gives a sys- 
tem with dimensionality less than 3; recent experiments 
in optical lattices have observed the properties of a one- 
dimensional Tonks gas in which bosons show fermionic 
properties. Several other experiments have realised con- 
ditions of two-dimensional confinement.i.e. where one 
degree of freedom is frozen out 0, 0, 0; however, the 
new physics in this regime remains to be explored: a 
two dimensional Bose gas in a homogenous potential does 
not undergo Bose-Einstein condensation (BEC), instead 
there is a Bcrcnzinskii-Kosterlitz-Thouless transition (a 
topological phase transition mediated by the spontaneous 
formation of vortex pairs), a system that is superfluid 
even though it docs not possess long-range order. This is 
counter the usual picture of superfluidity in three dimen- 
sions explained in terms of a macroscopic wavefunction 
describing the whole system. Early experiments on the 
KT transition were carried out with films of superfluid 
4 He 00 and more recent ones include the observation 
of quasi-condensates in thin layers of spin-polarized hy- 
drogen 0. 

In recent experiments, BECs created in conventional 
three-dimensional magnetic traps have been put into the 
quasi-2d (Q2D) regime through the addition of an opti- 
cal potential. In this limit the interaction energy, pro- 
portional to the chemical potential /i, is on the order 
of or smaller than the harmonic oscillator level spacing 
and along the tightly confined axial direction the char- 
acteristics of the condensate are those of an ideal gas . 
The crossover to the Q2D regime was first observed in 
and in by continuously removing atoms from a highly 
anisotropic trap to decrease the interaction energy. In 
the experiment described in the Q2D crossover is ob- 
served by gradually increasing the trap anisotropy from 
moderate to very large values whilst keeping the number 
of atoms fixed. 

The main aim of this paper is to use a simple theoret- 
ical model to examine some of the characteristic prop- 



erties of Q2D condensates, e.g. the chemical potential, 
release energy and quadrupole mode spectra are each cal- 
culated for a variety of various trap anisotropies which 
are achieved in the experiments. Note that a similar 
variational approach was taken in [8j which investigates 
the crossover to lower dimensions in general. We extend 
and expand upon this analysis by making a more general 
ansatz for the trial wavefunction (breaking the axial sym- 
metry) and deriving analytical expressions for chemical 
potential, total energy, release energy and radial conden- 
sate widths in the Q2D regime; a polynomial equation 
the solution of which gives the axial and radial conden- 
sate sizes in all regimes is also derived. It is used to 
calculate physical properties in all regimes and plot the 
crossover from the hydrodynamic to the Q2D limit for 
a number of physical quantities. The theoretical results 
are in very good agreement with the data of a recent ex- 
periment in 4]. We numerically calculate the frequency 
spectrum for the lowest quadrupole modes and also find 
an approximate analytic expression. 



The paper is structured as follows: In the second sec- 
tion we develop a formalism that allows us to calculate 
the ground state and the dynamics of a Q2D gas; start- 
ing from a Gaussian-parabolic trial wavefunction we use 
a variational method to obtain a set of differential equa- 
tions for the dynamics of the Q2D condensate and a poly- 
nomial equation for the ground state. We then examine 
the criteria for Q2D and calculate the chemical potential 
and energy per particle across the hydrodynamic-Q2D 
crossover. In section IV we calculate the size and re- 
lease energy of condensates in time-of-flight which are 
released from traps of varying anisotropy and make a 
quantitative comparison with recent experimental data 
0. The final section examines the change of collective 
excitation frequencies across the Q2D transition for the 
two quadrupole modes with zero angular momentum fol- 
lowed by the conclusions and a brief outlook. 
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II. THE HYBRID VARIATIONAL MODEL 



Condensates are usually trapped in harmonic poten- 
tials given by 



t r ^ 2 \ \ 2 2 



(1) 



where the (t) denote the trap anisotropics which can in 
general depend on time. A quasi-2D trap has X z \ x ,y 
For large anisotropics the condensate shape along the 
z-direction is very similar to the Gaussian profile of an 
ideal gas. However, along the weakly confined x- and 
y-axes the condensate has a parabolic shape characteris- 
tic of the hydrodynamic regime. The best description in 
terms of simple analytic functions is therefore to model 
the condensate wavefunction as a hybrid of parabola and 
Gaussian. Experiments on the condensate expansion in 
various regimes show the smooth crossover from hydrody- 
namic expansion to the characteristics of a quasi-2D gas, 
which essentially expands like an ideal non-interacting 
gas. To determine the dynamics of the quasi-2D conden- 
sate we use a variational method, as introduced in 
and define the trial wavefunction 



^ = A J l - 12 ~ ^e^e^+^y 2 ^ 2 ), (2) 



where the normalization constant A n is given by 



(3) 



The condensate width k(t) and phase Pi(t) parameters 
are functions of time and their time evolution completely 
describes that of the condensate. The condensate density 
profile is at all times restricted to a parabolic shape in 
the radial plane and a Gaussian shape along the highly 
compressed axial direction. The Lagrangian density for 
the nonlinear Schrodingcr equation is given by 



C = 



1 4 



-Vext M)H 2 + -5iV|v>| 



(4) 



with the nonlinearity parameter g = 47rfi. 2 a/m, where a 
is the scattering length, N is the number of atoms in the 
condensate and m is the atomic mass. After inserting 
the trial wave- function into Eq. Q the corresponding 
Lagrangian is found through integration L = J C d 3 x; 
the four terms of Eq.Q lead to 



V2gN 



(5) 



where we omitted the 'quantum pressure' term |Toj for 
the x and y directions (where this term is divergent due 
to the sharp boundaries of the condensate wavefunction 
in the hydrodynamic regime) but retained it for the z- 
direction where the condensate assumes the Gaussian 
shape of an ideal non-interacting gas (as the term pro- 
portional to f /I 2 ). The quantum pressure term is crucial 
in describing the dynamics. The total energy per particle 
E to t and the chemical potential fi are given by 



E 



kin ~f~-^pot 



-E u 



A* — £'kin+£ , pot + 2£ , int, (6) 



where -Ekin, E pot and E nrt are the kinetic, potential and 
interaction energy, given by the last three terms of the 
Lagrangian J3J), respectively. The Euler Lagrange equa- 
tions 

d_dL_dL d dL _ dL 

dt dk ~ dk ' dt ~ dfii ' ( ' 

yield the dynamic equations for the condensate widths k 
and phases /3j. We find for the widths 



k 



2h 



-Pik 



(8) 



After differentiating Eqs. (JSJ) once more w.r.t. time one 
can express the resulting second order equation in terms 
of the k alone: 
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(9) 

where Si Z — 1 for i — z and otherwise, ft is convenient 
to express the above equation in dimensionless quantities, 
so we introduce the dimensionless time r and widths di 
defined by 



A k 
di = — , 

ao 



(10) 



where ao = y/KJ {mujo) is the harmonic oscillator length. 
In terms of these quantities Eq. 1)10(1 can be rewritten as 

di = -XUt)di + - - f 1 - h lz ) + ±8 U , (11) 
dia x dyd z \ 3 J d% 



Li + L2 + £3 



R 2 1 2 3 2 l z 



La 



4Z? 



where the constant C v — 



±-2-N. To find the ground 

state of Eqs. fl 1(1 we have to set the left side equal to zero 
and solve the remaining coupled nonlinear equations: 



,2 y2 _ c v 
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(12) 
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This cannot be done analytically but it is straight forward 
to find a numerical solution. After some algebra and 
using various symmetries the three coupled equations can 
be reduced to one polynomial equation. Introducing new 
dimensionless units Di , defined as the condensate widths 
U normalized by the axial harmonic oscillator length a z , 
i.e. Di = U/a z , the polynomial equation can be written 
as 



8 1 

7 = - 

' 3 



r \ \ \ 1/2 

/ T" r j 



A 



3/2 



(13) 



where D z — j 2 . There is only one real and positive 
solution to this equation For the x- and y widths we find 



D, 



r a a 5/2 \ 1/4 



D, 



D = D — 

Uy — ^>x , 



(14) 



III. CRITERIA FOR QUASI-2D 

We shall now examine the case where the anisotropy 
becomes very large. A solution to Ea. (|13fl is then given 
by neglecting the first term on the r.h.s. and solving the 
remaining equation. We find that 7 2 = D z = 1 and thus 
the approximate solution is given by the axial harmonic 
oscillator length 



IzO — 



rriuj z 



(15) 



use the chemical potential shifted by an amount hui z /2 to 
calculate the radial width. Inserting solution lfT3|) for the 
axial width into Eq. 114(1 we obtain explicit expressions 
for the radial width 



l x0 




jy a ^yoKo 
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(17) 



and, after substituting into Eq. (fH)|) . for the chemical 
potential 



%to z 



^ N a ^ x0 ^ v ° 

n ao xi . 



(18) 



This expression shows that the chemical potential tends 
towards hui z /2, the harmonic oscillator ground state en- 
ergy, which is the energy per particle and also the chem- 
ical potential of the ideal non-interacting gas. The de- 
viation from this value is small for a quasi-2D gas, and 
interestingly, given by the same value as the correction 
to the axial condensate width of Ea. i|13|) . Note, how- 
ever, that this small deviation is vital for consistency. 
It is proportional to the square of the radial condensate 
width (see Eq. 116(1 . If it was zero the radial conden- 
sate width would also be zero which is neither possible 
nor self consistent. The expression for the Q2D chemical 
potential should be compared to that of a 3D hydrody- 
namic gas /i^d for which we obtain from |ll| after some 
rearrangements 



It is the minimum width the condensate shape can at- 
tain and it is also the solution for the width of a non- 
interacting gas. For this reason the gas along the z- 
direction is said to have the characteristics of an ideal 
non-interacting gas. We will see that this also applies 
to the expansion and the collective excitations of the gas 
which become identical to those of an ideal gas in the 
limit of large anisotropics. 

It is interesting to examine the range of validity of this 
approximation and find an estimate of the error. De- 
manding that the first term on the r.h.s. of Ea. l(13|) is 
much smaller than the second we see that the error or 
the deviation from the ideal gas solution scales with the 
ratio of NX xQ X y0 /X% 2 . The 2D -regime can be reached 
by either decreasing the number or increasing the axial 
frequency. 

Now we calculate the chemical potential from Eq.© 
and the terms of the Lagrangian JSJ and obtain after 
some algebra 



1 



„. 2;2 



x l xO> 



Tiuj z 



(16) 



where we used uj x I x0 = ujyl 2 Q , the expression for l z0 (Eq. 
115(1 and other symmetries of Eas. ((12|) . We find that the 
relation 1^ = y/2jl/mus 2 , i = x, y is similar to that of a 
hydrodynamic gas |ll| . only that for the quasi-2D gas we 



M3d 
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i/0 



(19) 



We observe that this expression tends towards zero for 
N — > and the power law is also different from the 
Q2D expression ((18|) . In previous papers [2| the con- 
dition fi3D < hu z was listed as a criterion for Q2D as 
H3D is on the order of the interaction energy at low den- 
sities. Similarly we can impose the conditon fi < Huj z 
on the expression of the Q2D chemical potential of Eq. 
((18(1 and we find for the maximum number of atoms to 
achieve 2D for a given trap geometry 



N < C 



9 \ ? ? ' 



(20) 



where the constant C — -y/32/225 for /13c < Tiio z and 
C = ^/ 7r/256 for fi < hw z . 

In Fig. we show plots of the chemical potential 
as a function of increasing radial trap frequency w r = 
lj x — Lu y , whereas the axial frequency remains constant 
at lo z /2t: = 2.2 kHz and the number of atoms is taken to 
be N — 8x 10 4 . For very small values of ui r the anisotropy 
A = X z /X x is very high and the chemical potential ap- 
proaches (i = huj z /2. Gradually increasing the radial 
trap frequency has the effect of reducing the anisotropy 
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contributions to the Lagrangian: 
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(21) 



Using Eq.© the total energy per particle is found from 
the Lagrangian (JSJ to be 



Etnt = 



2 



(22) 



For large anisotropics and small atom numbers this ex- 
pression tends towards Twj z /2, the harmonic oscillator 
ground state energy. This is what we expected, given 
that in the Q2D regime the gas is non-interacting and 
thus the energy per particle is equal to the chemical po- 
tential. The release energy, defined as the energy of the 
expanding cloud once the trap has been switched off, is 
given by the sum of the in-trap kinetic and interaction 
energies: 



E v 



E, 



l 



E u 



-J 1 



(23) 



which tends toward hu) z /4, equal to half the ground state 
energy, because the potential energy was lost when the 
atoms were released from the trap. We also find the 
atomic peak density ri pca k of the Q2D distribution, given 
by Eq. 0, from Eqs. P15I17| : 



FIG. 1: (a) The chemical potential plotted in units of hu z 
against increasing radial trap frequency: hybrid variational 
model (solid line), analytic approximation (dashed line), 
Gaussian variational model (dotted line); (b) the chemical 
potential plotted for a wider range of radial trap frequencies; 
the ideal gas limit is given by the dotted line, the hydrody- 
namic prediction by the dashed line. 



and increasing the chemical potential. An exact numer- 
ical solution using the variational hybrid model is given 
by the solid line. The dashed line denotes the analyti- 
cal approximation as given by Eq. i|18|) f° r a condensate 
well in the Q2D regime. There is very good agreement 
between the two for frequencies up to lj t ~ 20 Hz, where 
the anisotropy is on the order of A ~ 100 and /i « hui z . 
The dotted line is given by the results of a Gaussian vari- 
ational model as described in 0, where the trial wave- 
function in all spatial directions is given by a Gaussian. It 
is remarkable how similar the results are for both mod- 
els. In physical reality the actual wavefunction in the 
radial direction should be close to the inverted parabola 
we use in our 'hybrid model' (see Eq|5J. The plot in 
Fig. shows the chemical potential for a wider range 
of radial frequencies to r (solid line) on logarithmic scales; 
this shows how the values bend from the hydrodynamic 
asymptote (dashed line) to the ideal gas value of %u> z /2 
(dotted line), with the bending point given by values of 
u> r /2ir f« 20 Hz and [i ~ huo z . From Eqs. (|5lt)ll6|l we find 
the following useful relations between the various energy 



"peak = NA 2 n = - = 



(24) 



IV. CONDENSATE EXPANSION AND 
RELEASE ENERGY 



It is interesting to examine the expansion of the con- 
densate in time-of-flight (TOF) in the 2D-regime. As 
a Q2D gas resembles an ideal gas in the axial direc- 
tion its expansion is described by the Schrodinger equa- 
tion and can be calculated analytically. Looking at Eqs. 
©, noting that X z (t) = for expansion and neglect- 
ing the second term describing the interaction we obtain 
h = h 2 /(m 2 ll). The solution to this second order differ- 
ential equation for l z (0) — I z q and l z = is given by 




= hoVl +L0 2 z t 2 , 



(25) 



where Eq. (|15J) has been used. In the following we will re- 
fer to this simple expression as the 'ideal gas expansion'. 
Fig.[5]shows the axial and radial condensate widths, nor- 
malized by their initial values (in the trap), during TOF 
for a trap with N = 8 x 10 4 atoms, uj z /2tt = 2.2 kHz 
and w r /27r = 7 Hz, which are typical parameters for the 
experiment in 4J. The prediction of the hybrid model, 
which is indistinguishable from that of the Gaussian vari- 
ational model, is given by the solid line. For the axial di- 
rection, shown in (a), it is very close to the expansion of 
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FIG. 2: The condensate axial (a) and radial (b) widths, nor- 
malized by their initial values, plotted in TOF. The solid line 
is based upon the hybrid model, the dotted upon the expan- 
sion of a non-interacting wavepacket and the dashed line upon 
the hydrodynamic theory. 



the ideal gas given by Eq. (|2*5|l and plotted as the dotted 
line. The faster expansion of a hydrodynamic gas 01 
is given by the dashed line. These curves demonstrate 
clearly that for the axial expansion the system behaves 
like an ideal gas. This contrasts the radial expansion 
of Fig. 03, where the hybrid and Gaussian variational 
models (solid line) agree nearly perfectly with the hydro- 
dynamic theory (dashed line) but are far removed from 
the ideal gas expansion (dotted line); this demonstrates 
the hydrodynamic character of the radial expansion. 

In Fig. |3 we plotted the condensate expansion as a 
function of time; however, in the experiment in Q the 
expansion times are kept constant at 15 ms and the radial 
trap frequency u> r is varied to explore the various regimes 
and the crossover to Q2D. The measurements in £j were 
done for two different axial trap frequencies and atom 
numbers. The results of the measurements are shown in 
Fig-EK where the open and filled circles are the data for 
traps in which the atoms initially have axial oscillation 
frequencies of iv z /2ir — 1990 Hz and 960 Hz, respectively. 
To obtain a theoretical comparison we propagate the hy- 
brid model Eqs. i|ll|) repeatedly for varying initial con- 
ditions obtained from Eq. lfT3|) . The hybrid model pre- 
dictions for the two optical traps are given by the solid 
lines and we find good agreement to the experimental 
data for both traps. The horizontal dashed lines indicate 
the expansion of the ideal gas, given by Eq. (|25[> and 
the dotted lines the expansion of the hydrodynamic gas. 
The 'ideal gas' and hydrodynamic models yield straight 
lines on the log- log plot. In contrast, the curve for the hy- 
brid model follows the hydrodynamic asymptote down to 
uj r w 20 Hz, corresponding to [i rj hu> Zl where it bends 
and follows the 'ideal gas' line towards zero radial fre- 
quency. This transition from the hydrodynamic to the 
'ideal gas' asymptote gives conclusive evidence of the gas 
entering the Q2D regime. There is a noticeable devia- 
tion for u> r well below the cross-over frequency because 



in the experiment imperfections broke the radial sym- 
metry and added a residual potential along the y-axis. 
This was taken into account in the theoretical calcula- 
tions. There is good agreement with the experimental 
data. The overall frequencies in x- and y-direction are 

given by uj x — uj r and uj v = + (Wy es ) 2 . The residual 

frequencies ujy CS were measured in Q and are given by 
u T y es /2Tr = 26 Hz for the trap with uj z /2tt = 1990 Hz and 
10^/27: = 12 Hz for the trap with u z /2n = 960 Hz. 

It is also possible to calculate the release energy E vc \ of 
the condensate from the experimental measurements of 
its axial width as has been done in |2( . The release energy 
is easily calculated after a long time-of-flight. In this case 
the potential term in Eq. (JSJd) is zero and the interaction 
term negligible so that we are left with the kinetic energy 
alone and find E Ic i = E^ in . The terms in x and y can 
be neglected as most energy is in the previously tightly 
confined z-direction and we obtain from the Lagrangian 
(j3J) an expression for the release energy 

h 2 ( , 9 1 \ m -o 1 ft 2 1 , % 
m \ Ali I 4 4 to 1% 



rel 



^kin 



where we used relation (jHJ) for the last step. After 
an initial acceleration, when released from the trap, 
the condensate moves with constant velocity and after 
long enough time-of-flight (TOF) we can approximate 
l z « l z /t, where t is the TOF. The release energy is then 
written as 

mil , 
Er^^f, (27) 

where we neglected the quantum pressure proportional 
to 1/7 2 as it tends to zero for long time-of-flights. In 
a Q2D gas the total energy is given by the sum of the 
potential and the kinetic energy in the axial direction, 
both of which contribute by an equal amount to the har- 
monic oscillator ground state. As the potential energy is 
lost when the trap is switched off the remaining kinetic 
energy constitutes half the ground state energy so that 
_E rc i = Tiu z /A. The release energies for the same set of 
data as in (a) are plotted in Fig. |3Jd, Ea. H27|) was used to 
determine the release energy from the condensate width. 
The theoretical predictions of the hybrid model are given 
by the solid lines, the dashed lines indicate the value for 
the ground state and the first excited state, respectively. 
In 0] the expression Ef el = mL 2 /(14i 2 ) is used to find 
the release energy, where L z is the width of a fitted in- 
verted parabola. This formula can be derived by making 
a fully parabolic ansatz to the trial wavefunction - an 
approach that is valid in the hydrodynamic limit but is 
not really applicable in the Q2D limit, where one should 
rather use Eq. H27|) as otherwise the release energy is 
underestimated by about 20%. 

There is excellent agreement between the theoretical 
predictions and the experimental data and we find that 
the release energy does indeed go towards Tiui z /A indicat- 
ing that the Q2D limit was reached in the experiment 
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FIG. 3: (a) The axial expansion of the condensate in 
TOF for two different trap geometries and atom num- 
bers. Solid lines indicate the theoretical predictions, dashed 
lines indicate the ideal gas limit and dotted lines the hy- 
drodynamic limit. The data are taken for traps with 
initial oscillation frequencies (before release) of u> z /2n = 
1990 Hz (open circles), 960 Hz (filled circles) Q. The atom 
numbers are 8 x 10 4 and 1.1 x 10 5 for the upper and lower 
curves, respectively, (b) release energies for the same data 
points; the dashed lines indicate the energies of the ground 
and the first excited state. 



V. COLLECTIVE EXCITATIONS IN Q2D 



Another, as yet unexplored, method to observe the 
transition to Q2D is to probe the collective excitation 
spectrum. An ansatz for the trial wavefunction of the 
type (J2J allows for the description of the three quadrupo- 
lar modes llj. In an axially symmetric trap they are 
given by the m = 2 mode and the m = low- and high- 
lying modes, where m denotes the angular momentum 
quantum number. In order to calculate the mode fre- 
quencies we linearise the dynamic equations of the hybrid 
model (|llf) around the ground state. Making the ansatz 



(28) 



inserting it into Eas. (|ll|) . expanding up to first order in 
e, and using Eas. (|12fl to simplify and combine certain 



terms, we obtain after some algebra 
( 



3C P 




yO 
Cp/3 



yO 
C P /3 



d lo d y« d lo 



\ d to d yO d io d x® d yO d io ~ d lrO d yO d ^o d to I 



e z 



(29) 

Calculating the eigenvalues and eigenfrequencies of the 
above matrix one finds the collective excitation frequen- 
cies and modes. This can be easily done numerically. For 
simplicity that lends itself to an easy analytical treatment 
we assume cylindrical symmetry X x o = X v q = A. r o, d x o = 
dyo ~ 
of two 



d r Q and reduce the set of three equations to a set 



Dd z0 



tO 



(30) 



where D — C p /(d^. dl ). We can diagonalize the above 
matrix and find the eigenvalues which yields for the eigen- 
frequencies uj 2 : 



~ 2 — 2A^ — -Ddro + 2 

Wn 6 




-2X 2 z0 + -Dd r0 + 2 



2 

-D 2 d 2 

3 20 



(31) 



Inserting expression (|15I17|I for the ground state widths 
into the equation above we obtain an analytic expres- 
sion for the collective excitation frequencies in the Q2D 
regime: 



9X 2 1 C pKo , o 



A 



zl) 



± 



\ 




2 y C P X%X^ . 



The two frequencies given by Eq. l|31|l describe the high- 
and low-lying m = modes of the collective excita- 
tion spectrum. The high-lying m — mode is an in- 
phase compressional mode along all directions (breath- 
ing mode). The low- lying m = corresponds to a radial 
oscillation of the width which is out of phase with an 
oscillation along the trap axis. The third mode (not de- 
scribed by Eq. 13 1 1> in an axially symmetric trap is the 
m = 2 mode. It corresponds to a quadrupole type ex- 
citation in the radial plane and its frequency is given 
uo2 = \/2k-V, irrespective of the axial frequency. Fig. 0] 
shows the change of the collective excitation frequencies 
of the two m = modes for increasing u> r . We find 
that the high- lying mode frequency (u>+) changes from a 
value of lj + = 2 iu z for radial frequencies close to zero 
to a value of d + = \/3 uj z for large radial frequencies. 
The latter value has been determined from the hydrody- 
namic model |ll| . given by the dashed line, in the limit 
of uj r — ► 0. The hardly distinguishable predictions of the 
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FIG. 4: The m = high- lying mode frequency and 
the m — low-lying mode frequency (u)-) plotted against 
the radial trap frequency u r /2n in figures (a) and (b), re- 
spectively; hybrid variational model (solid line), Gaussian 
variational model (dotted line) and hydrodynamic prediction 
(dashed line). 

hybrid- and Gaussian variational models are given by the 
solid and dotted lines, respectively. For very small u r , in 
the Q2D regime, the mode frequencies approach the ideal 
or non-interacting gas limit where the mode frequency is 
twice the trap frequency. In the ideal gas limit the radial 
oscillation goes towards zero and we obtain a pure axial 
oscillation for the high-lying mode, which can be found 
from the eigenvectors of (|29|) . 

Something similar happens for the low-lying mode fre- 
quency w_ which changes from the hydrodynamic limit 
of cj_ = y/10/3 0J r to oj- = 2co r for decreasing uj r . An 
analysis of the eigenvector (corresponding to this mode) 
of matrix l|29[) shows that in the Q2D regime the axial 
component of the oscillation is increasingly suppressed 
and goes towards zero in the limit of infinitely small lo t . 



Surprisingly, in this limit the oscillation in the radial di- 
rection at the mode frequency w_ = 2uj r seems not at 
all affected by the hydrodynamic character and strong 
interactions of the gas in the radial plane. It oscillates at 
the same frequency as an ideal collisionless gas, although 
it is far from being collisionless. This feature has been 
pointed out by Pitaevskii and Rosch and depends 
on the special symmetry in the two-dimensional regime. 



VI. CONCLUSIONS AND OUTLOOK 

We studied the properties of Q2D condensates using a 
hybrid variational model based on a Gaussian-parabolic 
trial wavefunction. The chemical potential and ground 
state energy were calculated for a wide range of parame- 
ters from the hydrodynamic to the Q2D regime. We find 
that the chemical potential approaches the harmonic os- 
cillator ground state and derive an analytical expression 
for its values in the Q2D regime. The condensate size in 
time-of-flight and its release energy are plotted for traps 
with varying anisotropy. We find good agreement to re- 
cent experimental results We also calculate the ex- 
citation spectrum of the quadrupole modes of the Q2D 
gas and find a gradual change from the hydrodynamic 
values to values equal to twice the trap frequencies in 
the Q2D regime, as predicted in All results from 

the hybrid model are compared to those derived from a 
Gaussian variational model |9(, a hydrodynamic model 
[pH] and the ideal gas model. In future work we want to 
explore the effect Q2D has on vortex structure, dynamics 
and decay. 
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